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1.  Typical 
bubbler  system 
arrangement. 


Site 

conditions 


completely,  creating  an  open-water  area  above  the  bubbler, 
or  it  can  reduce  its  thickness,  so  that  ships  can  pass  through  it 
or  so  the  ice  breaks  up  earlier  and  more  easily. 

Bubblers  transfer  heat  to  the  ice  cover  through  the  circula¬ 
tion  of  water  caused  by  a  rising  plume  of  air  bubbles  released 
from  orifices  well  beneath  the  surface  (Fig.  1).  The  water 
from  the  lower  levels,  usually  warmer  than  that  near  the  ice 
cover,  brings  heat  to  the  ice  and  causes  melting.  Under  favor¬ 
able  conditions,  the  energy  going  into  melting  can  be  many 
times  greater  than  the  energy  used  to  operate  the  system.  It  is 
this  multiplication  effect  that  makes  air  bubblers  so  advanta¬ 
geous  over  any  direct  use  of  energy  to  melt  ice. 

However,  an  air  bubbler  system  can  operate  successfully 
only  when  the  characteristics  of  the  site  are  favorable.  Several 
factors  have  to  be  known  before  any  bubbler  installation  can 
be  planned. 

It  should  be  obvious  that  the  physical  dimensions  of  the 
water  body  need  to  be  known.  A  fairly  good  idea  of  the  depth 
and  bottom  topography  (plus  the  type  of  bottom  materials, 
for  anchoring  purposes)  is  necessary.  The  depth  and  the  hori¬ 
zontal  dimensions,  when  combined  with  the  temperature  of 
the  water,  give  an  estimate  of  the  heat  that  can  be  drawn  on 
to  melt  the  ice. 

Water  temperature  is  of  prime  importance.  It  can’t  just  be 
assumed  that  the  temperature  at  the  bottom  is  4°C  (39  °F), 
the  maximum-density  temperature.  While  theoretically  this 
would  be  expected,  actual  temperature  measurements  of 
water  bodies  in  winter  (which  are  unfortunately  very  few) 
have  shown  that  in  most  water  bodies  where  bubblers  might 
be  installed,  actual  bottom  water  temperatures  are  much  less, 
as  low  as  0.1°  to  0.2  °C  (32.2°  to  32.4  °F)  or  even  0.01  °C 
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(32.02  °F).  Field  investigations  indicate  that  only  in  very  deep 
lakes  or  very  small,  sheltered  lakes  do  bottom  temperatures 
approach  4°C  (39  °F).  And  in  a  water  body  with  any  signifi¬ 
cant  amount  of  flow  or  circulation,  such  as  a  river  or  a  bay 
open  to  a  larger  water  body,  the  temperature  may  be  at  0°C 
(32  °F)  throughout  the  depth  as  a  result  of  turbulent  mixing, 
although  temperatures  as  high  as  0.5°  to  1.0°C  (32.9°  to 
33.8  °F)  have  been  found  by  Wortley  in  some  bays  of  the 
Great  Lakes  (G.  Ashton,  personal  communication,  1982). 

Natural  movement  of  the  water,  through  either  flow  or  cir¬ 
culation,  is  important  to  evaluate,  not  only  because  of  the 
chance  of  having  no  available  reserve  of  heat,  but  also  be¬ 
cause  the  flow  will  affect  the  water  circulation  created  by  the 
bubbler.  The  expected  location  of  melting  can  be  shifted 
horizontally,  with  the  result  that  the  actual  amount  of  warm¬ 
er  water  moved  up  to  the  ice  cover  becomes  diluted  and  less 
effective.  Naturally,  flow  is  expected  in  rivers.  But  flow  may 
be  present  in  lakes  or  harbors  too,  due  to  the  inflow  of  tribu¬ 
tary  streams,  wind-induced  circulation,  or  currents  from  a 
large  water  body  connected  to  a  bay  or  harbor. 

To  evaluate  the  feasibility  of  a  bubbler  system,  some 
knowledge  of  the  natural  ice  cover  characteristics  is  needed. 
For  example,  the  normal  ice  thickness  should  be  known,  so 
that  it  can  be  predicted  whether  the  bubbler-induced  melting 
will  create  open  water  or  simply  thin  the  ice  to  some  equilib¬ 
rium  thickness.  In  the  case  of  open  water,  the  expected  size  of 
the  opening  is  related  to  the  normal  ice  thickness.  Also,  it  is 
critically  important  to  know  if  the  ice  cover  will  be  stationary 
or  moving.  If  it  is  moving,  the  melting  will  be  spread  over  a 
large  area  of  ice  but  there  may  be  so  slight  a  reduction  in  ice 
thickness  as  to  be  of  no  benefit.  Expected  dates  of  ice  cover 
formation  and  break-up  should  also  be  determined. 

Finally,  the  information  about  site  conditions  should  in¬ 
clude  climatic  data.  The  air  temperature  at  the  site  will  influ¬ 
ence  the  heat  transfer  through  the  ice  cover.  The  amount  of 
snow  on  the  ice  cover  also  governs  the  heat  transfer  rate,  be¬ 
cause  its  insulating  effects  will  keep  the  top  surface  of  the  ice 
warmer  than  the  air  if  the  air  temperature  is  below  0°C 
(32  °F).  The  best  indicator  of  temperatures  affecting  an  ice 
cover  is  the  freezing  index,  the  total  number  of  degree-days  of 
freezing.  The  more  detailed  the  freezing  index  information 
the  better,  e.g.  average  cumulative  totals  of  degree-days  dur¬ 
ing  the  winter,  or  consecutive  daily  records.  Useful  additional 
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operation 


climatic  information  includes  wind  data,  cloud  cover  data, 
precipitation  records,  etc.,  although  such  data  are  not  readily 
factored  into  the  currently  available  mathematical  analysis  of 
air  bubbler  operation. 

Mechanical  Details.  The  size  and  type  of  air  compressor 
used  in  a  bubbler  installation  depend  on  many  factors.  The 
size  and  depth  of  the  bubbler  will  dictate  the  air  require¬ 
ments.  Whether  the  installation  is  permanent  or  not,  and 
where  it  is  located,  will  help  to  determine  the  type  of  com¬ 
pressor  used.  Compressors  commonly  fall  in  the  range  of 
0.0094  m’/s  (20  ft’/min)  to  0.71  m’/s  (1500  ft’/min)  in  deliv¬ 
ery  rates,  at  pressures  up  to  3.4  x  10s  Pa  (50  psig).  It  is  some¬ 
times  determined  that  an  air  dryer  is  needed  in  the  system,  to 
prevent  ice  formation  at  the  orifices.  In  order  for  this  deci¬ 
sion  to  be  made,  the  expected  relative  humidity  of  the  air  at 
the  site,  plus  the  pressure  drop  at  the  orifices,  will  have  to  be 
considered. 

Piping  from  the  compressor  is  made  up  of  the  supply  line 
and  the  diffuser  line.  Standard  design  approaches  govern  the 
choice  of  the  supply  line,  balancing  economy  and  perform¬ 
ance.  Naturally,  it  is  helpful  to  keep  the  supply  line  short  to 
minimize  frictional  losses,  but  the  most  convenient  location 
for  the  compressor  station  will  probably  be  the  main  factor 
setting  the  supply  line’s  length.  So  its  diameter  should  be 
chosen  to  permit  pressures  high  enough  at  the  end  of  the  sup¬ 
ply  line  to  ensure  air  discharge  at  the  most  distant  orifice  of 
the  diffuser  line.  A  rapid-emptying  fitting  is  often  supplied  at 
the  far  end  of  the  diffuser  line,  to  permit  rapid  removal  of 
water  at  start-up. 

The  diffuser  line  may  present  a  design  challenge  in  some  in¬ 
stallations.  There  is  a  step-wise  variation  in  air  flow  and  vel¬ 
ocity  as  each  orifice  is  passed,  and  there  is  a  gradual  variation 
in  pressure  along  the  length  of  the  line.  As  suggested  above, 
air  must  reach  the  orifice  farthest  downstream,  and  pressure 
has  to  be  higher  than  the  water  pressure  at  the  level  of  the 
line.  In  addition,  a  fairly  constant  air  discharge  rate  per  unit 
length  of  line  is  desired  from  one  end  to  the  other.  For  instal¬ 
lations  that  are  large  or  that  have  nonuniform  orifice  spacing 
or  diffuser  pipe  sizes,  a  detailed  analysis  of  the  diffuser  is 
needed  (see  Vigander  and  others,  1970).  But  for  most  installa¬ 
tions  of  modest  size,  the  variations  in  air  discharge  along  the 
pipe  are  negligible.  Berlamont  and  others  (1973)  and  Camp 
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and  others  (1968)  give  techniques  that  evaluate  variations 
with  fixed  orifice  spacing  and  pipe  size. 

As  a  general  rule,  orifice  spacing  is  chosen  to  be  half  the 
depth  of  submergence  of  the  diffuser.  This  spacing  gives  rea¬ 
sonably  good  results  approximating  a  line-source  of  air,  as 
opposed  to  separated  point-sources.  Orifice  size  depends  on 
the  desired  air  discharge  rate,  and  is  determined  by  the  con¬ 
ventional  orifice  discharge  equation: 


sQ a  =  Qo  =  Cd^J-V2Ap/ea 


where  s 

Qa 

Qo 

Cd 

d 

Ap 


Pa 


orifice  spacing 

air  discharge  rate  per  unit  diffuser  length 
air  discharge  rate  from  one  orifice 
orifice  discharge  coefficient 
orifice  diameter 

difference  between  air  pressure  within  the  dif¬ 
fuser  and  water  pressure  outside  the  diffuser 
mass  density  of  the  air  within  the  diffuser. 


Plume  Mechanics  and  Melting.  Considering  the  bubbler  as 
a  line  source  of  air,  the  plume  it  creates  is  two-dimensional. 
Figure  2  shows  the  velocity  distribution  of  the  plume  (a  Gaus¬ 
sian  distribution),  which  increases  linearly  with  height  above 
the  vicinity  of  the  diffuser  line.  Because  of  flow  conditions 
near  the  orifices  and  the  initial  momentum  of  the  air  jet,  the 


2.  Definition 
sketch  for  a 
two-dimensional 
bubbler  plume. 
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linear  boundaries  of  the  plume  have  a  theoretical  vertex 
about  0.8  m  (2.6  ft)  below  the  diffuser  (denoted  by  x0  in  Fig. 
2).  The  expansion  of  the  plume  results  from  drawing  in 
steadily  increasing  amounts  of  surrounding  water  (entrain¬ 
ment),  all  driven  by  the  momentum  of  the  buoyant  rising 
bubbles.  It  is  this  entrainment  that  delivers  heat  from  the  pre¬ 
sumably  warmer  water  at  depth  to  the  underside  of  the  ice, 
and  causes  melting.  The  water  flow  in  the  plume  per  unit 
length  (£?w)  at  the  surface  can  be  described  in  terms  of  Qa,  air 
discharge  per  unit  length,  and  H,  the  depth  of  the  diffuser: 

<2w  =  5.5302  S(//  +  0.8)°-5[loge(l+  -j^j)]05 

which  can  be  approximated  by: 

Qw  =  2.28 

This  is  in  meter-second  (SI)  units,  and  is  for  sea  level,  though 
the  effect  of  other  elevations  is  minor. 

As  shown  in  Figure  2,  the  plume  is  2b  wide  when  it  reaches 
the  surface.  The  measure  of  b  is  chosen  as  the  standard  devia¬ 
tion  of  the  Gaussian  velocity  distribution  of  the  plume  at  the 
surface.  So  the  width  of  the  plume  in  meters  can  be  found  by: 

2b  =  0.364(//  +  0.8)gJ15. 

After  the  plume  strikes  the  underside  of  the  ice,  or  reaches 
an  open  water  surface,  the  flow  spreads  horizontally  outward 
from  the  area  of  impingement.  It  is  not  important  for  this  dis¬ 
cussion  to  examine  this  horizontal  flow,  other  than  to  say 
that  its  velocity  drops  off  very  rapidly,  and  the  principal  melt¬ 
ing  effects  of  the  bubbler  occur  in  and  near  the  impingement 
region  of  width  2b  as  defined  above. 

What  is  important  is  what  happens  to  the  horizontal  flow 
after  that,  i.e.  the  induced  circulation  within  the  water  body. 
For  the  bubbler  to  continue  to  melt  ice  over  a  long  period,  it 
must  draw  on  the  thermal  reserve  not  just  in  the  vicinity  of 
the  bubbler,  but  over  a  broad  area.  Its  ability  to  do  this 
depends  on  the  circulation  induced  in  the  water  body,  a  sub¬ 
ject  that  is  not  well  understood. 

Two  possibilities  have  been  envisioned,  a  “filling”  scheme 
and  an  “overturning”  scheme.  In  the  filling  model,  the  pre- 
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sumably  less  dense  water  resulting  from  bubbler  operation 
accumulates  at  the  surface,  gradually  displacing  the  denser 
fluid  downward.  In  the  overturning  model,  the  plume  sets  up 
circulation  cells,  with  the  plume’s  flow  plunging  back  to  the 
level  of  the  diffuser  line  at  some  distance  to  the  side  of  the 
bubbler.  Work  by  Baines  and  Turner  (1969)  indicates  that  cir¬ 
culation  takes  place  in  water  bodies  whose  depth  is  greater 
than  about  one-fourth  their  width,  and  that  the  filling  model 
is  more  likely  to  be  followed  in  shallower  or  broader  water 
bodies. 

Recent  studies  (G.  Ashton,  personal  communication,  1982) 
suggest  that  circulation  and  complete  mixing  take  place  with¬ 
in  a  horizontal  distance  of  5  to  7  times  the  depth  for  any  bub¬ 
bler  depth,  regardless  of  the  width  and  depth  characteristics 
of  the  water  body.  Outside  this  circulation  zone,  it  appears 
that  exchange  of  water  still  occurs  via  layered  flow  toward 
and  away  from  the  circulation  zone.  Thus  a  combination  fill¬ 
ing  and  overturning  scheme  forms  a  third  possible  model.  In 
a  strict  filling  model,  the  thermal  reserve  would  be  used  by 
the  bubbler  at  an  initial  high  rate  which  decreases  with  time  as 
the  thermal  reserve  is  used  up.  In  a  strict  overturning  model 
the  circulation  cells  would  prevent  heat  from  being  drawn  in¬ 
to  the  bubbler  from  beyond  the  first  cell.  Thus  the  actual 
available  thermal  reserve  would  be  much  less  and  would  be 
depleted  much  earlier.  In  a  combination  model,  the  thermal 
reserve  consumption  would  lie  between  these  two  extremes. 

How  much  heat  is  delivered  to  the  surface  by  the  bubbler? 
This  of  course  depends  on  the  water  flow  delivered  to  the  sur¬ 
face,  Qw,  and  on  the  temperature  profile  in  the  water  from 
the  surface  down  to  the  depth  of  the  diffuser.  Although  not 
strictly  accurate,  it  is  an  acceptable  approximation  to  take  the 
temperature  of  the  plume  as  it  reaches  the  area  of  impinge¬ 
ment  as  the  depth-weighted  average  of  the  temperature  pro¬ 
file  between  the  surface  and  the  diffuser  line,  7~avg.  Then  the 
heat  delivered  to  the  surface  is 

&  =  Gwl^avg”  ^Cw^p 

where  G  =  heat  delivered  per  unit  time  per  unit  length  of  the 
two-dimensional  bubbler 
=  plume  flow  at  the  surface  (described  earlier) 

Ts  =  temperature  at  the  surface  (0  °C  in  the  case  of  an 
ice  cover) 


on 
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ew  =  mass  density  of  water 

cp  =  water’s  specific  heat  capacity. 

It  must  be  remembered  that  the  temperature  profile,  and 
hence  T^g,  change  with  time  as  the  bubbler  continues  to  op¬ 
erate. 

As  most  of  the  heat  G  remains  with  the  water  in  circula¬ 
tion,  we  need  to  know  how  much  is  actually  transferred  to  the 
ice  cover.  By  using  an  analogous  analysis  for  a  slot  jet  of  air 
impinging  on  a  plate,  Ashton  developed  relationships  to  yield 
the  heat  transfer  coefficient  within  the  impingement  region. 
In  terms  of  Qa  and  H,  the  heat  transfer  coefficient  in  W  m'2 
°C'‘  (SI  units)  is 

hb  =  17,220[!oge(//+  10.3)  -2.332]°‘3l(//->-0.8)-0'69Qa'16 
which  can  be  approximated  by 


/iu  =  6888  H 


-  0.384  n0. 16 
Via 


This  coefficient  is  considered  constant  over  the  width  2b,  and 
drops  off  very  rapidly  to  either  side  of  the  bubbler.  Thus  the 
actual  heat  transferred  to  the  ice  cover  (in  W  m*2)  in  the  im¬ 
pingement  region  is 

db  =  ^b^^avg  “  TO 


db  ~  ^b'^avg 

since  Ts  =  0°C  for  an  ice  cover. 

What  does  this  mean  in  terms  of  thinning  an  ice  cover  or 
maintaining  an  ice  cover  at  an  equilibrium  thickness?  The  an¬ 
swer  lies  in  the  basic  heat  balance  equation  at  the  ice/water 
interface, 

db~di  =  Q\Hdr}/d() 

where  qx  =  heat  conduction  upward  through  the  ice  (and 
through  a  snow  cover  if  present)  to  the  atmos¬ 
phere 

Ci  =  mass  density  of  the  ice 
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X  =  latent  heat  of  melting 
dt)/dt  =  melting  rate  of  the  ice  of  thickness  r?. 

If  qb  is  greater  than  qb  then  the  ice  cover  will  become  thinner, 
while  if  qb  =  qt,  it  will  maintain  an  equilibrium  thickness. 

The  calculation  of  qt  is  straightforward,  and  won’t  be  dis¬ 
cussed  here.  It  should  be  noted,  however,  that  a  snow  cover 
will  reduce  qx  significantly,  and  that  snow  cover  thicknesses 
and  the  conductivities  of  ice  and  snow  will  need  to  be  known. 

Two  examples  may  serve  to  illustrate  the  application  of  the 
foregoing  information.  First,  suppose  a  bubbler  is  capable  of 
transferring  qb  =  200  W  m"2  to  the  ice  cover,  and  that  the 
mean  daily  air  temperature  (an  acceptable  substitute  for  the 
temperature  of  the  top  surface  of  a  snow-free  ice  cover)  is 
-25°C.  The  equilibrium  ice  thickness  jjc  would  be  found  by 
the  relation 


db  = 


HT5-Tm) 

% 


where  k{  =  2.24  W  m"'  deg'1,  the  conductivity  of  pure  ice;  Ts 
=  -25°C,  the  top  ice  surface  temperature;  and  =  0°C, 
the  bottom  ice  surface  temperature.  Then  i)e  is  found  to  be 
0.28  m.  (If  instead  qb  was  400  W  m"2  or  100  W  m‘2,  then  ije 
would  be  0.14  m  or  0.56  m,  respectively.) 

The  second  example  considers  a  thinning  ice  cover.  Sup¬ 
pose  Qa  =  0.001  m2  s'1  and  H  =  10  m,  such  that  hb  =  979  J 
nr2  s'1  deg  ',  and  Tavg  =  0.1  °C.  Then  qb  =  97.9  W  m'2. 
Now,  for  simplicity,  assume  that  the  top  surface  ice  tempera¬ 
ture  is  0°C,  so  that  q,  =  0  W  nr2,  and  the  entire  amount  of  qb 
goes  to  melting.  Then  qb  =  Qt\dr)/dt.  Taking  gj  =  916  kg 
nr’  and  X  =  3.34 xlO5  J  kg"',  then  di\/dt  becomes  3.12xl0'7 
m  s"'  or  2.76  cm  per  day  of  melting,  a  bit  more  than  an  inch  a 
day. 

Rationale  of  Operation.  Any  prospective  bubbler  installa¬ 
tion  should  first  be  subject  to  a  calculation  of  the  maximum 
theoretical  thermal  reserve  available  for  melting  ice.  Then 
this  should  be  translated  into  the  maximum  theoretical  vol¬ 
ume  of  ice  melting.  What  this  will  show,  in  many  cases,  is 
that  the  amount  of  ice  suppression  that  can  be  achieved  is  not 
as  extensive  as  expected,  or  in  other  words,  that  the  thermal 
reserve  is  a  scarcer  commodity  than  first  thought.  This  scarce 
thermal  reserve  must  be  used  cautiously  and  prudently. 
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Two  operating  schemes  can  be  chosen,  singly  or  combined, 
to  wisely  exploit  the  reserve.  First,  instead  of  seeking  the  crea¬ 
tion  of  open  water,  most  system  objectives  can  be  met  just  as 
well  by  achieving  a  thinned  ice  cover.  This  greatly  reduces 
heat  wasted  to  the  atmosphere,  and  helps  to  save  the  thermal 
reserve  so  that  it  can  be  drawn  upon  over  a  longer  period  of 
time.  Similarly,  intermittent  operation  of  the  bubbler  system, 
rather  than  continuous  operation,  will  extend  the  usefulness 
of  the  thermal  reserve.  Qualitatively  matching  periods  of 
bubbler  operation  to  times  when  air  temperatures  produce 
the  greatest  potential  ice  growth  would  be  particularly  benefi¬ 
cial,  but  even  an  arbitrary  on-off  schedule  is  helpful. 

The  performance  of  a  bubbler  system  is  dependent  on  the 
air  flow  rate  Qa  and  the  depth  of  submergence  H ,  both  of 
which  the  designer  can  choose,  plus  the  water  temperature 
and  thermal  reserve,  which  he  can’t.  To  illustrate  the  ex¬ 
pected  values  of  some  of  the  important  parameters,  Table  1 
gives  values  of  £)w,  the  water  flow  per  unit  length  reaching 
the  impingement  region,  in  m2  s'1  (ft2  min"');  2b,  the  width  of 
the  impingement  region,  in  m  (ft);  and  hb,  the  heat  transfer 
coefficient,  in  W  m"2  °C~'  (Btu  ft"2  s'1  °F"');  all  for  three  typi¬ 
cal  values  of  air  flow  Qa  and  five  typical  values  of  depth  of 
submergence  H.  The  same  information  (except  heat  transfer 
coefficient)  is  portrayed  in  Figure  3. 

Tabic  1.  Bubbler  performance  parameters  for  three  air  flow  rates  (g,)  and 
five  submergence  depths  (//). 


tm) 

tfo 

tm'  s') 

(ft*  min  ') 

tm) 

tft) 

tWm  '  °CV 

f Bluff's''  °F  ) 

Q,  = 
2 

0.01  ml  9 

7 

"  <6.5  fl*  min  •> 

0.39  250 

0.51 

1.7 

2.37x10' 

3.76  x  10  ' 

6 

20 

0.98 

630 

1.2 

3.9 

1.73x10' 

2.74x10" 

10 

33 

1.5 

970 

2.0 

6.6 

1.42  x  10' 

2.25x10" 

20 

66 

2.6 

1700 

3.8 

12 

1.04  x  10‘ 

1.65x10  ’ 

40 

130 

4.4 

2800 

7.4 

24 

7.37  x  10' 

1.17x10" 

H 

.  ™ 

& 

0.001  m* 

7 

s'1  <0.65  ft'  win'1) 

0.12  78 

0.36 

1.2 

1 C4  x  10* 

2.60  x  10" 

6 

20 

0.31 

200 

0.88 

2.9 

l.20x  10' 

l.90x  10 

10 

33 

0.47 

300 

1.4 

4.6 

9.80  x  10-’ 

1.55x10" 

20 

66 

0.83 

540 

2.7 

8.9 

7.19x  10' 

1.14x10 

40 

130 

1.4 

900 

5.3 

17 

5.l0x  10* 

8.08  x  10 

0.  ■ 
2 

0.0001  m 

7 

1  s  '  <0.065  ft1  min  ') 

0.039  25 

0.26 

0.85 

1 .14  x  10' 

1.81  x  10 

6 

20 

0.098 

63 

0.62 

2.0 

8.27  x  10' 

1.31  x  10 

10 

33 

0.15 

97 

0.99 

3.2 

6.78  x  10- 

1.07  x  10 

20 

66 

0.26 

170 

1.9 

6.2 

4.97  x  10* 

7.88  x  10  ? 

40 

130 

0.44 

280 

3.7 

12 

3.53  x  10 

5.60  x  10 

•  • 


>  * 


Melting  Ice  With  Air  Bubblers 


11 


0Q  (m2s_l) 


Ashton,  G.D.  (1974)  Air  bubbler  systems  to  suppress  ice. 
USA  Cold  Regions  Research  and  Engineering  Laboratory, 
Special  Report  210,  37  pp. 

Ashton,  G.D.  (1978)  Numerical  simulation  of  air  bubbler  sys¬ 
tems.  Canadian  Journal  of  Civil  Engineering,  5:  231-238. 
Ashton,  G.D.  (1979)  Point  source  bubbler  systems  to  sup¬ 
press  ice.  Cold  Regions  Science  and  Technology,  1:  93-100. 
Baines,  W.D.  and  J.S.  Turner  (1969)  Turbulent  buoyant  con¬ 
vection  from  a  source  in  a  confined  region.  Journal  of  Fluid 
Mechanics,  37(Part  I):  51-80. 

Berlamont,  J.  and  A.  Van  der  Beken  (1973)  Solutions  for 
lateral  outflow  in  perforated  conduits.  Journal  of  the  Hy¬ 
draulics  Division,  American  Society  of  Civil  Engineers, 
99(HY9):  Proceedings  Paper  10004,  pp.  1531-1549. 

Camp,  T.R.  and  S.D.  Graber  (1968)  Dispersion  conduits. 
Journal  of  the  Sanitary  Division,  ASCE,  94(SA1):  Proceed¬ 
ings  Paper  5781,  pp.  31-39. 

Vigander,  S„  R.A.  Elder  and  N.A.  Brooks  (1970)  Internal 
hydraulics  of  thermal  discharge  diffusers.  Journal  of  the  Hy¬ 
draulics  Division,  ASCE,  96(HY2):  Proceedings  Paper  7085, 
pp.  509-528. 


3.  Water  flow  (Q^ 
m2  s' ',  dashed 
lines)  and  impinge¬ 
ment  width  (2b,  m, 
solid  lines )  as  func¬ 
tions  of  air  dis¬ 
charge  and  bubbler 
depth. 


References 


FILMED 


